Metallic mean-field stripes, incommensurability and chemical potential in cuprates 
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We perform a systematic slave-boson mean-field analysis of the three-band model for cuprates with 
first-principle parameters. Contrary to widespread believe based on earlier mean-field computations 
low doping stripes have a linear density close to 1/2 added hole per lattice constant. We find 
a dimensional crossover from ID to 2D at doping ~ 0.1 followed by a breaking of particle-hole 
symmetry around doping 1/8 as doping increases. Our results explain in a simple way the behavior 
of the chemical potential, the magnetic incommensurability, and transport experiments as a function 
of doping. Bond centered and site-centered stripes become degenerate for small overdoping. 

PACS numbers: 71.28.+d,71.10.-w,74.72.-h,71.45.1r 
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It is now a well established fact that doped holes in 
cuprates tend to self-organize in antiferromagnetic (AF) 
domain walls ^, |^, ^ ^ . These quasi one-dimensional 
(ID) structures called stripes where predicted by mean- 
field theories 1^, Q inspired by the problem of solitons in 
conducting polymers^. 

In the superconducting systems observed stripes are 
parallel to the Cu-0 bond (hereafter called vertical 
stripes) and for low doping have a linear filling fraction 
J/ « 1/2 of added holes per lattice constant (a = 1). Un- 
less a one-dimensional (ID) instability opens a gap at the 
Fermi level^, half-filled {f = 1/2) stripes are expected 
to be metallic. On the contrary mean-field theories pre- 
dicted ly =: I insulating stripes which led to an early 
rejection of stripes [p^. 

Here we perform a systematic analysis of slave-boson 
(SB) mean- field jlll solutions for a three-band Hubbard 
model as appropriate for the cuprates. Contrary to the 
early mean-field analysis we find that for low doping 
X in a material like La2-a;Sra;Cu04 (LSCO) the most 
favorable stripes have indeed v ^ 1/2 (and are verti- 
cal) reconciling the mean- field picture with experiment. 
Small doping introduces more stripes at fixed u reduc- 
ing the stripe distance d and explaining the observed be- 
havior of the magnetic incommensurability e « x with 
6 = l/(2d)|,|, 1 1,§. 

We find low doping stripes to be centered on oxy- 
gen, bridging two Cu sites, hereafter named bond cen- 
tered (BC) stripes after the terminology of one-band 
models |l|, |l| (see Fig.|). 

Significant interstripe overlap sets in at a; ~ 0.1 trigger- 
ing a crossover from ID behavior at low doping to 2D be- 
havior at large doping. This does not affect the stability 
of stripes up to x = 1/8 where a change of regime occurs. 
Further doping occurs at a fixed stripe distance d = 4 and 

starts to grow. This explains the change of behavior 
observed in both the incommensurability|^, |[ D and 
the chemical potential[^ |l^ around doping x — 1/8. 
Furthermore for a; < 1/8 the chemical potential is at 
the center of an approximately particle-hole symmetric 
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FIG. 1: BC vertical stripe with d = 7 and = 1/2 (a; = 
1/14). Sites with (without) arrows are Cu (O). The length of 
the arrows is proportional to the local magnetization on Cu 
(negligibly on O). The radius of the circle represent the excess 
charge respect to the AF state (0.09 holes for the largest one). 
The y axis is taken on the stripe direction (vertical). The 
picture repeats periodically in x and y directions 



band (Fig. ^ whereas for doping x > 1/8 particle-hole 
symmetry is broken in good agreement with the picture 
deduced from transport experiments p^, p7| . Finally at 
X w 0.21 BC stripes and site-centered (SC) stripes be- 
come degenerate suggesting a regime of soft lateral stripe 
fiuctuations possibly relevant for superconductivity. This 
quasidegeneracy was found in more sophisticated treat- 
ments of (less realistic) one-band mode ls [p^ , 13 however 
the ordering of the states was reversed [|l3|. 



One-band models can provide qualitative understand- 
ing but may miss subtle but important details specific 
for the cuprates. Therefore we use a three-band Hub- 
bard model with a parameter set deduced from a first- 
principle LDA computation for LSCO|l^. We view our 
computation as the last stage of a first principle compu- 
tation and thus without free parameters. The largest 
repulsion C/d|l8) is treated within a SB approach [pT| 
fully equivalent to the Gutzwiller approximation and the 
other interactions are treated within the Hartree-Fock 
(HF) approximation (an equivalent formalism has been 
used in Ref. [|l^ for electron-doped cuprates). We solve 
the mean-field equations both in real and in momen- 
tum space. The real space computations are completely 
unrestrictedly except for the spins which are assumed 
to be along the z axis. 

Typically we minimize the energy in large clusters with 
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FIG. 2; Electron mean-field bands for the stripe of Fig. |l| 
measured from the chemical potential. Left (right) panel is 
in the direction perpendicular (parallel) to the stripe. We 
also plot the insulating charge transfer gap at momentum 
k = (0,7r) measured from the same reference energy (dots). 

Ns stripes {Ns — 2,4) of a maximum length L = 30. 
Depending on boundary conditions, size and filling, so- 
lutions with and without charge order in the stripe di- 
rection are obtained. In general due to the ID character 
of stripes one expects that a ID instability will develop 
in large systems at low doping rendering the system in- 
sulating even for v < I. Indeed if superconductivity is 
suppressed by pulsed magnetic fields the system is insu- 
lating below a critical doping and a crossover tempera- 
ture r* |2l|| . In this work for simplicity we concentrate 
on the metallic state and restrict to solutions without 
charge order in the stripe direction. Since T* is quite 
small compared with the typical electronic scales we can 
set the temperature to zero in our computations. 

Fig. 1^ depicts the typical charge and spin distribution 
for a BC vertical stripe. Added charge is accumulated 
in the antiphase boundary between the AF regions and 
mainly resides on oxygen sites. 

In Fig. ^ we show the band structure for the electrons 
in an extended zone scheme. Two bands appear well in- 
side the charge transfer gap. The upper one at ~ 0.9 eV 
is more Cu like and quite flat (even in the stripe direc- 
tion) whereas the lower one has substantial dispersion, 
is mainly O like, and crosses the Fermi level. Counting 
both ingap bands the ly — 1/2 system corresponds to 3/4 
hole (1/4 electron) band filling. The gap between the 
bands is related to the magnetic character of the stripe. 
(It is absent in SC solutions with a non-magnetic core.) 
The electronic structure close to the chemical potential 
is well represented by a half-filled cosine like band, here- 
after referred to as the active band. 

The dispersion of the active band perpendicular to the 
stripe is quite flat indicating quasi ID behavior. Small 
oscillations appear in this direction for d < 5 which can 
be identified as a crossover interstripe distance from ID 
to 2D behavior and corresponds to a crossover doping 
X = 0.1 for V = 0.5. 




0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 
Stripe filling fraction v 







p 




</=87 


A 




(b) 


1 hole U \ ^ 

>K n \ 5 










« \ 

\ 0.1 5 
1^^' ^ 0.1 
0.05 


'/ ^ 


Theory - 
-xCuO,^, O 
^Sr,Cii04 V 







La, ,,_^Ndy 


jSr^CiiOj >S 




0.1 X 


0,2 



0.05 0.1 1/80.15 0.2 0.25 0.3 
Hole doping x 

FIG. 3: Ch as a function of v (a) and x (b) for vertical BC 
stripes. Sizes are 4dx 12 for d < 6 and 2dx 12 for d > 6. In (a) 
we also show the result for 4d x 30 (open squares) and setting 
tpp — shifted down by 0.89 eV (crosses). In (b) we also show 
the result for site-centered stripes (labeled SC, open symbols) 
for d = 3, 4 and one self-trapped hole or "electronic polaron" 
a,t X — 1/36. The inset in (b) reports the incommensurability 
as obtained from the present calculation (line) compared with 
experimental data from Ref. [^. 

The band structure in the stripe direction (right panel) 
is quite similar in all BC stripes solutions regardless of 
d which instead determines the periodicity in the per- 
pendicular direction (left panel). Comparison with pho- 
toemission data requires matrix elements considerations 
and will be presented elsewhere. Here we notice that the 
Fermi surface crossing at vertical momentum ky w 7r/4 
has been recently observed for d = 4 stripes j2^ and the 
flat portion of the active band in the direction perpen- 
dicular to the stripe correlates well with the flat bands 
observed in many cuprates around k « (tt, 0)[p3[. 

To evaluate the stability we compute the energy per 
added hole for Nh holes added to the system with N 
Cu02 units, eh = [E{N + Nh) - E{N)]/Nh. Here 
E{N + Nh) [E{N)] is the total energy of the doped [un- 
doped AF] solution. In Fig. ^ (a) we show eh as a func- 
tion of the flUing fraction v = Nh/{NsL) for BC verti- 
cal stripes. Each curve corresponds to a fixed d. The 
curves have a sharp minimum at v = 0.5. This is in con- 
trast to early one-band ||] and three-band mean-field 
computations for which v = 1 is the most favorable 
hole filling. We can reproduce this behavior by setting 
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artificially tpp = in our parameter set (crosses). The 
difference in behavior can be understood by analyzing 
the band structure. In the case tpp — the active band 
becomes as narrow as the 0.9eV band (Fig. ||) and a gap 
develops between the active band and the band imme- 
diately bellow instead of overlapping as in Fig. ||. This 
gap produces the cusp singularity a,t v — I [Fig. || (a)] 
via the discontinuity in dE/dN which makes v — 1 very 
stable. As tpp is increased the gap closes, the cusp singu- 
larity disappears, and a minimum develops in the v < 1 
part of the curve. A simple model explains the forma- 
tion of this minimum. Consider a flat density of states 
(DOS) of width W to estimate the kinetic energy of the 
active band and neglect interactions. In this approxima- 
tion e/i = cq/v + Wv/2 where eo is the energy cost per 
unit length to create a domain wall in the undoped sys- 
tem (eo ~ J/2 with J the superexchange constant). The 
optimum filling is given by Vo-pt — min(-\/2eo/W, 1). In- 
creasing tpp increases W without substantially affecting 
eo and produces a local minimum for v < 1. The stabil- 
ity close to = 1/2 is further enhanced because the DOS 
of a ID band is not flat but has a minimum at v = 1/2. 
For the same reason a gap produces a cusp in the energy, 
a minimum of the DOS enhances the local curvature of 
eh vs. u and tends to shift the optimum filling closer to 
V = 1/2. 

We have checked that for the most common parameter 
sets which contain tpp[|l8| results are similar. However, as 
soon as tpp is set to zero, as in Ref. half-filled stripes 
become unstable in favor of filled ones. 

The tpp = result is in accord with HF computations 
in the simplest one-band model. Including a second- 
neighbor hopping {t') we have found that partially filled 
stripes can he stabilized also in one-band models by a 
similar mechanism, which will be reported elsewhere. 

Partially filled stripes where also found in the one-band 
Hubbard model with i' = in a much stronger correla- 
tion regime and going beyond a static mean-field In 
our case, because the active band is broad, correlations 
are less important and we expect that a static mean-field 
provides a good starting point. 

Vertical stripes should be compared with other pos- 
sible ground states. We find that within the present 
SB approach they are lower in energy than polaron 
solutions Q [see Fig. ^ (b)] and diagonal stripes solu- 
tions (not shown). Very low nonsuperconducting dopings 
{x < 0.05) were not explored in detail since we believe 
that a careful consideration of other effects is required in 
this case. Especially long-range order perpendicular to 
the planes favors loop configurations of stripes [jl^ and 
long-range Coulomb effects!^, ^ are expected to be- 
come important. 

If one uses the HF approximation instead of SB the 
minimum also occurs at ^ 1/2 for BC vertical stripes. 
However a polaron lattice is the ground state in HF and 
diagonal stripes are lower in energy than vertical ones. 



Thus the overall success of the present mean-field com- 
putation with respect to the earlier ones|^ is due to both 
a more accurate mean-field approximation SB instead of 
HF) and a more accurate parameter set. 

Further stabilization of the half-filled stripes can occur 
if due to many-body effects a gap or pseudogap tends to 
open at the commensurate filling i' — 1/2^. Although 
we do not need this effect to explain v 1/2 it is quite 
possible that this produces a fine tuning for T < T*. In 
this regard it is interesting to remark that the "V" shape 
form of the curves in Fig. |3| (a) is due to the gap produced 
by the discretization of the levels in a finite system. For 
larger systems [see d = 4 curve in Fig. ^ (a)] the cusp 
becomes rounded, however the minimum is still close to 
v = 1/2. One should be aware that the curves with the 
cusp are for already quite long stripes {L = 12) and it is 
not clear whether stripes in real materials will be much 
longer. Thus finite size data could turn out to be more 
realistic than infinite size one. 

For d > 4 all curves coincide close to ly — 0.5 in 
Fig. H (a) whereas for d = 3,4 the curves are shifted 
up. This feature is due to the width of the domain wall 
of 4 5 lattice sites (see Fig. |l|) which forces stripes to 
overlap leading to an increase of e^- 

In Fig. H (b) we show et for various values of the stripe 
separation d and as a function of doping x assuming each 
dopant introduces one hole. Since x = v/d curves with 
larger d appear at lower concentration. The locus of the 
minimum of eh as a function of doping in Fig. ^(b) is ex- 
pected to form a continuous curve in the thermodynamic 
limit by combining different d solutions in new solutions 
with larger periodicity. Up to x w 1/8 the stripe fill- 
ing is fixed at 1/ « 0.5 and consequently the density of 
stripes increases with doping. This explains the behav- 
ior of the incommensurability e = x/(2h') x as seen in 
neutron scattering experiments in this doping range (in- 
set of Fig. |)|l|, |, |, |, |. For a; > 1/8 the right branch 
of the d = 4 solution is more stable than the v « 0.5 and 
d — 3 solution due to the stripe overlap effect. There- 
fore the incommensurability remains locked at e « 1/8 in 
good agreement with the change of behavior in e observed 
aroundx^ l/8g|||. 

As doping increases further, BC stripes become de- 
generate with SC ones at a; « 0.21. This suggests that 
as doping increases lateral fluctuations of the stripe will 
become soft and possibly mediate pairing between holes. 

For doping x > 0.225 we find the d = 3 stripe 
(e = 1/6 « 0.17) to become the lowest energy solution 
with an initial filling fraction v — 0.675. Experimen- 
tally the situation is not clear, e = 1/6 has been re- 
ported for YBasCuaOe+^ll (YBCO) but not for LSCO 
where e remains in the e = 1/8 line up to a; = 0.25 (in- 
set of Fig. I). In La2-x-yNdySr:rCu04 (LNSCO), where 
stripes are pinned by the low-temperature tetragonal lat- 
tice distortion, the incommensurability is substantially 
increased beyond e = 1/8 at doping a; = 0.2 but without 
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reaching e = 

It is possible that lateral stripe fluctuations become so 
strong at this doping range that an effectively isotropic 
state is reached |27|. Another possibility is that the 
d = 3 phase is skipped due to phase separation among 
the d = A phase and the overdoped Fermi liquid which for 
our parameters becomes the lowest energy solution close 
0.4 (see Ref. 
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for a related scenario). Fur- 
ther theoretical and experimental work should be done 
to clarify this point. Especially being the system charged 
a careful analysis of phase separation is needed p^. 

It has been recently emphasized that there is a close 
connection between the doping dependent incommensu- 
rability as discussed above and the chemical potential 
in cuprates. In fact the chemical potential for the elec- 
trons can be related to Ch via: ^jl — — (e/i -I- x^^). From 
Fig. ^(b) one can then deduce that /i is approximately 
constant for x < 0.1 and decreases for x > 0.1 in qualita- 
tive agreement with the observed behavior 1 14, |l^. The 



rate of change of /i with doping, being a high derivative 
of the energy, is very sensitive to finite size effects and, 
moreover, few experimental points are available in this 
doping range in order to allow for a precise comparison. 
A rough estimate indicates that the theoretical rate of 
change of with doping for a; > 1/8 is approximately a 
factor of 2 larger than the experimental onejl^, [Tst . This 
may be due to phase separation among the d ~ A stripe 
solution and the paramagnetic overdoped Fermi liquid as 
mention above. 

We finally turn to the discussion of our calculation in 
light of recent measurements of the Hall coefficient Rjj in 
LNSCO|l§ and YBC0(l3. In the former compound 
displays an abrupt decrease below the charge-stripe or- 
dering temperature Tq and for concentrations x < 1/8. It 
has convincingly been argued that quasi ID transport is 
not enough to explain this anomaly and instead reflects a 
remarkable cancellation due to particle-hole symmetry in 
the stripe state[^^ |2^. A partial suppression of R// (and 
simultaneously the thermopower S) below some tempera- 
ture Tmax has also been observed in YBCO up to oxygen 
contents corresponding to doping a; ~ 1/8 For dis- 
tant stripes we indeed observe that the chemical potential 
crosses an approximately particle-hole symmetric band 
(see Fig. Interstripe hopping becomes significant for 
X > 0.1 however it does not break particle- hole sym- 
metry. As discussed above the chemical potential shifts 
from the center of the band for x > 1/8, thus breaking 
particle-hole symmetry and consequently both Rh and S 
start to grow in modulus. Unfortunately the magnitude 
and even the sign of the transport coefficients cannot be 
predicted by a knowledge of the band structure alone. 
For example Rh for our Fermi surface (which turns out 
to be open) is determined by the unknown anisotropy of 
the scattering path length 

To conclude we have shown that, contrary to 
widespread believe, an appropriate mean-field theory of 



stripe phases in cuprates is not in contradiction with 
experiments. Our theory explains in a simple way the 
change of behavior of several experimental quantities (in- 
commensurability, chemical potential, transport coeffi- 
cients) around the "magic" doping x — 1/8. Obviously 
the present stripes are static whereas those in real materi- 
als are dynamic unless pinned by a structural distortion. 
Many properties however will be largely independent of 
the dynamical or static character of the stripes. The 
energy for example is determined by short range correla- 
tions and is insensitive to the long-distant behavior. This 
is analogous to the problem of spin waves in a ID anti- 
ferromagnet where the spin wave approximation gives a 
reasonable value for the energy although the system does 
not have true long-range order Since the properties 
we have discussed are mainly determined by the ener- 
getics we believe that the present mean-field theory is a 
good starting point for an understanding of the electronic 
structure of hole doped cuprates. 
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